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Abstract. We consider a quadratic matrix boundary value problem with equations and 
boundary conditions dependent on a spectral parameter. We study an inverse problem that 
consists in recovering the differential pencil by the so-called Weyl matrix. We obtain asymptotic 
formulas for the solutions of the considered matrix equation. Using the ideas of the method of 
spectral mappings, we prove the uniqueness theorem for this inverse problem. 
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O. 

CN . 1. Introduction and main results 



In this paper, we consider the boundary value problem L = L(£, U, V) for the equation 

£Y :=Y" + (p 2 ■I + 2ipQ 1 (x)+Q (x))Y = 0, x £ (0, 7r), (1) 



p i ■ with the boundary conditions 

^ ■ U{Y) := Y' (0) + {iph x + h )Y(0) = 0, 

V(Y) := Y'{tt) + (ipH! + H Q )Y(n) = 0. 



(2) 



Here Y(x) = [yk{x)] k= Y^ is a column vector, p is the spectral parameter, / is the m x m unit 
matrix, Q s (x) = [Q s ,jk( x )]j,k=T^n are m x m matrices with entries Q s ,jk(x) G H^ffO, 7r], s = 0, 1 
h s = [h s ,jk]j,k=T^> H s = [H s ,jk\j,k=T^i> where h s>jk , H s>jk are complex numbers. 

We assume that det(7 ± hi) ^ and det(J ± Hi) ^ 0. This condition excludes problems of 
f-^. \ Regge type (see from consideration, as they require a separate investigation. 

Differential equations with nonlinear dependence on the spectral parameter, or with so-called 
"energy-dependent" coefficients, frequently appear in mathematics and applications (see J2HH] 
^ . and references therein). In particular, inverse problems for such equations arise in investiga- 
tion of mathematical models describing collisions of relativistic spinless particles [7] or proper 
vibrations of a string with viscous drag [8]. 

In this paper, we investigate the inverse problem for the pencil L, which consists in recovering 
coefficients of the boundary value problem (II]), ([2]) by its spectral characteristics. In the scalar 
case (m = 1) inverse problems for quadratic pencils were studied in works jHHIl]- 

In the particular case when Qi(x) = 0, equation ([I]) becomes the matrix Sturm-Liouville 
equation. In recent years, significant progress has been made in the inverse problems theory 
for this equation. Constructive algorithms for solution of inverse problems were suggested, and 
characterization of spectral data was given (see [ToTTTS] ). 

Now we plan to use ideas, developed for the matrix Sturm-Liouville equation, for problems 
with nonlinear dependence on the spectral parameter. In the present paper, we consider a 
general situation, without any conditions of selfadjointness on the coefficients and with the 
spectral parameter in the boundary conditions. We study the inverse problem for the pencil L 
by so-called Weyl matrix and prove the uniqueness theorem for the solution of this problem. 
For our investigation, we develop ideas of the method of spectral mappings [20PT] . This method 
also can be used to obtain a constructive procedure for the solution of this inverse problem, 
but this question requires separate investigation. 
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One of the main difficulties in problems for differential pencils is related to asymptotic 
behavior of solutions. The main terms of asymptotic representations depend on the coefficients 
of the pencil Qi, hi, Hi. Derivation of asymptotic formulas is nontrivial even in the scalar case, 
and in the matrix case it is more complicated, because of additional difficulties connected with 
noncommutativity of matrix multiplication. We obtain the asymptotics for the fundamental 
system of solutions for equation flTJ). They are required to prove uniqueness, but they also can 
be considered as a separate result. 

Now proceed to the formulation of the main result. Let <3>(x,p) = [$jfc(x, p)]j t k=T^m be the 
matrix solution of equation (JTJ) satisfying the conditions U($) = I, V(&) = 0. We call $(x,p) 
the Weyl solution for L. Put M(p) := $(0,p). The matrix M(p) = [Mjk(p)]j tk= i^ is called 
the Weyl matrix for L. The notion of the Weyl matrix is a generalization of the notion of the 
Weyl function (m-function) for the scalar case (see J201E2]) and the notion of the Weyl matrix 
for the matrix Sturm-Liouville operator (see [T81IT9] ). 

The inverse problem is stated as follows. 

Inverse Problem 1. Given a Weyl matrix M(p), construct the coefficients of the pencil 

L. 

Remark. One can prove in a standard way that the boundary problem L has a countable 
set of eigenvalues {p n }- The entries of M(p) are meromorphic in p and their poles coincide 
with {p n }- As in the scalar case (see [llj), the following representation is valid 



where m n are multiplicities of the corresponding eigenvalues p n , and M nv are some matrix 
coefficients. 

Following [TT1JT2"] . we call the collection {p n , M nu } the spectral data of the pencil L. By virtue 
of ([3]), the spectral data determine the Weyl matrix uniquely. Therefore, Inverse Problem 1 is 
equivalent to 

Inverse Problem 2. Given spectral data {p n , M nu }, construct the coefficients of the pencil 

L. 

In this paper, we restrict ourselves to Inverse Problem 1. 

Along with L we consider a pencil L of the same form but with other coefficients Q s (x), h s , 
H s . We agree that if a symbol 7 denotes an object related to L then 7 denotes the corresponding 
object related to L. 

We now state the uniqueness theorem for Inverse Problem 1. 

Theorem 1. If M(p) = M(p), then L — L. Hence the Weyl matrix determines the coefficients 
of the pencil ([T]), (TSJ) uniquely. 

In order to prove Theorem 1, we need asymptotics for the solutions of equation ([T]). They 
are obtained in Section 2. In Section 3, we provide the proof of the uniqueness theorem. 

2. Asympotic behavior of the solutions 

Let C(x,p) = [Cjk(x, p)]j,k=T~m an d S(x,p) = [Sjk(x, p)]jk=T^m be the matrix solutions of 
equation (TTJ) under the initial conditions 



The main goal of this section is to obtain the asymptotics of S(x,p) and C(x,p) as \p\ — > 00. 




(3) 



C(0,p) = S'(0,p)=I, 



C"(0,p) = S(0,p) = 0. 
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Let the matrix functions P + (x), P-(x), P*(x) and Pl(x) be the solutions of the Cauchy 
problems 

P! t (x) = ±Q 1 (x)P ± (x) 1 P±(0)=/, () 
P£{x) = ±P£(x)Q l (x), P|(0) = I. 1 ' 

Remark. In the scalar case (m = 1), one has 

P±0) =P±{x) =exp|±^ Qi(t)cft 

Lemma 1. T/ie following relations hold 

P + (x)P*(x) = Pl(x)P + (x) = I, P_(x)P*(x) = = /, 

for all x G [0, 7r]. 
Proof. Using (j3J), we get 

(P:(.t)P + (x))' = (P!(x))'P + (x) + P!(x)P|(x) = -P*(x)Q!(x)P+(x) + Pl(x)Q 1 (x)P + (x) = 0. 

Hence P*(x)P + (a;) does not depend on x, so P*(x)P + (x) = P*(0)P+(0) = I. The other 
relations can be proved similarly. □ 

Theorem 2. For x e [0, 7r], |p| — > oo ; z/ = 0,l, i/ie following relations hold 

C {u) (x,p) = ( ^exp(tpx)P4x) + ( ^^exp(-ipx)P + (x) + Oip"- 1 exp(\r\x)), (5) 
SM(x,p) = ^ 1 exp(zpx)P_(x) + 1 exp(-zpx)P + (x) +0(p^ 2 exp(|r|x)), 

where r := Imp. 

In order to prove the theorem, we develop the ideas of [23] with necessary modifications 
caused by the matrix case. 

Proof. 1. First we derive Volterra integral equations for S(x,p) and C(x,p). 
One can easily check that the matrix functions 

C (x,p) = ^exp(ipx)P_(x) + - exp(-ipx)P + (x), 

S (x,p) = J—exp(ipx)P-(x) - — exp(-ipx)P + (x) 
Zip Zip 

form a fundamental system of solutions for the differential equation 

Y" + Q\{x){ip ■ I - Q 1 {x))~ 1 Y' + (p 2 • I + 2ipQ 1 (x) - Q\{x))Y = 0. 

Rewrite equation ([1]) in the form 

Y" + Qi(x)(zp • / - Q 1 (x))- 1 Y' + (p 2 • I + 2i / oQ 1 (a;) - Q 2 (x))F = P(x,p, Y"), (6) 
F(x,p, Y) := Q[(x)(ip ■ I - Qi{x))- l Y' - {Q\{x) + Q (x))Y. (7) 
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Apply the method of variation of parameters to this equation. Every solution of 
represented in the form 

Y(x, p) = C (x, p)A(x, p) + S (x, p)B(x, p), 
where coefficient matrices A(x,p) and B(x,p) can be found from the system 



can be 



C (x 


P) 


S (x 


P) ' 




' A 1 









C' Q (x 


P) 


Sq(x 


P) _ 




B' 


[x,p) _ 




_F{x,p,Y) _ 



and the initial conditions 

A(0,p) = Y(0,p), B(0,p) = {I-Q l (0)/{ip))- 1 Y'(0,p). 
In order to solve system (181) , we introduce the matrix functions 

Co(x,p) = ^exp(ipx)P*(x) + ^exp(-ipx)P*(x), 
Sq(x,p) = — exp(ipx)Pl(x) - — exp(-ipx)P*(z). 

Using Lemma [H one can easily show that 



(8) 



(9) 



Therefore 



" Co 


So 


-1 




-S* ' 




ip{ip ■ 


I - 


Qi)- 1 o 


_C 






-ci 


Co . 









ip(ip ■ I - QiY 1 






A(x,p) = - 


-ipSo(x, 


p){ip-I - 


Qi 


»)- 1 F(x,p,F), 



B'(x, p) = ipC*(x, p)(ip ■ I - Q l {x))- l F{x, p, Y). 



and 



Y{x, p) = C {x, p)A(0, p) + S {x, p)P(0, p)+ 

+ ip [ X (S (x, p)C*(t, p) - C (x, p)S*(t, p))(ip ■ I - Q 1 (t))- 1 F(t, p, Y) dt. (10) 
Jo 

Simple calculations show that 

S Q (x,p)C*(t,p)-C {x,p)S*(t,p) = ^{exp(ip(x-t))P_(a;)P;(t)-ex P (-ip(x-t))P + (a;)P*(t)}. 

Substituting Y = C(x,p) and Y = S(x,p) into ([9]) and ffTOl) . we arrive at the following integral 
equations 

C(x,p) = C (x,p) + X - J X {ex P (tp(x - t))P_(x)P;(t)- 

- exp(-ip(z - t))P+(x)Pl(t)}(ip ■ I - Qtit^F^p, C) dt, (11) 

gi(0) ^ 1 + ~ J\eMip^ - t))P-(x)Pl(t)- 



S(x,p) = S (x,p) I 



ip 



exp(-ip(x - t))P + (x)Pl(t)}{ip ■ I - gi(0) _1 P(t,p, S) dt, 
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where F(x,p,Y) is defined in ([7]). 

2. Then we continue to work with C(x,p). The function S(x,p) can be treated similarly. 
Differentiating (fTT|) with respect to x and using fll]), we get 

C'(x,p) = C' (x,p) + l -{ip ■ I - Q 1 (x)) J X {ex P ( ip (x - t))P_(ar)P;(t)+ 

+ exp(-zp(x - i))P+(x)P*(t)}(zp • / - Qtf^Ffap, C) dt, (12) 

Denote 

p y (p) := max ||C^(ar, p) exp(— |r|x) ||, 1/ = 0, 1. 

XS[0,7r] 

We agree to denote by the same symbol K different positive constants not depending on 
x and p, and we mean by |.|j the following matrix norm: ||v4|| = max J= j^ J2T=i \ a ij\i A = 

[ a jk]j,k=l,m- 

Since 

\\F(t,p,C)\\ < K\\Q[(t)\\^- + \\Ql(t) + Q (t)\\M, 

\p\ 

we get from ffTTl) and (TL2l) 



M <k(i + ^ + ^§-\, „ lW < k ( ipi + , M + ^ 

whence we obtain p u (p) < K\p\ u or 

C {u) (x,p) = 0(|p| v exp(|r|a;)). 
Substituting this into ffTTj) and (fT2|) . we arrive at ([5]). □ 

3. Proof of the uniqueness theorem 

Let <p(x,p) = [<pjk(x, p)]j } k=T~m an< i = [4>jk(x, p)]j,k=T~m be the matrix solutions of 

equation ([!]) under the initial conditions 

<p(0, p) = ^(tt, p) = /, f%) = Vf» = 0. 

It is easy to check that 

p) = S(x, p) + <p(x, p)M(p), (13) 
$(x,p)=^(a;,p)(f/(^))- 1 . (14) 
We expand ip(x,p) by the fundamental system of solutions C(x,p) and S(x,p) 

cp(x, p) = C(x, p) - S(x, p)(iph 1 + h ) 

and, using Theorem [2J we get 

tpM (x,p) = ^ exp (tpx)P_ (x) (I-h) + exp (-tpx)P + (x) (I + h\)+ 

+ C^p^ 1 exp(|r|a;)), x G [0, 7r], v = 0, 1, |p| ->■ oo. 

In order to obtain the asymptotics for ip(x,p), we can apply the substitution x — )■ 7r — x, 
hi — > — /io - ► — -f^o- Then we need the analogs of P±(x) and P±(x), which we denote by 
P±{x) and P±*(x) and define as the solutions of the Cauchy problems 

P^(x) = tQi(x)P^x), P£(7r)=/, 
The following lemma establishes connections between these asymptotic coefficients. 
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Lemma 2. For x G [0, ir], the following relations hold 

P + (x)Pl(0) = P'(x), P-(x)Pl(0) = Pl(x), 
P'*(0)P;(x) = P'*(x), Pl*(0)Pl(x) = P+{x). 

The proof, based on using the corresponding Cauchy problems, is trivial. 
One can obtain the following asymptotic formula 

^\x,p) = t^expMvr - x))Pl(x)(I + H x ) + exp(-ip(7T - x))P' + (x)(I - H x )+ 

+ 0(p^ 1 exp(|r|(7r-x))), xE [0,vr], u = 0, 1, |p| -> oo, 

and then using (flij) . one can derive asymptotic formulas for Q(x, p) and M(p). 

Hereafter the asymptotics in the angles 9^ := {p G C: 5 < ±argp <n — 5},0<5<n 
will be required, so we formulate the following result. 

Lemma 3. Suppose x G (0, it), v — 0, 1, |p| — > oo; then 
(i) for p G 6^~, we have 

^\x,p) = -t^exp(- W )P + (x)(/ + / il ) + 0(p^ 1 exp(|r|x)), 

P) = exp(-zp(7r - x))P;(x)(J - H 1 ) + exp(|r|(7r - x))), 

= -(-ip)"- 1 exp(ipx)P;(x)(P;(0))- 1 (/ + /i 1 )- 1 + 0(p^ 2 exp(-|r|x)), 
M(p) = (^p)- 1 (/ + / il )- 1 + 0(p- 2 ). 

fwj /or p G 6^, we /jai>e 

tp("\x,p) = ^^exp(zpx)P_(x)(/-/i 1 ) + 0(p^ 1 exp(|r|x)), 

</> H (^P) = ^^exp(zp(7r - x))Pl(x){I + H x ) + O^" 1 exp(|r|(7r - z))), 

$H(x,p) = -{ipy- 1 exp{-tpx)P'{x){P'{0))-\l - h)- 1 + 0{p"~ 2 exp{-\r\x)), 
M( P ) = -(tp)-\l-h 1 )- 1 + 0(p- 2 ). 

By definition, put 

tZ := Z" + Z(p 2 ■ / + 2ipQ 1 (x) + Q {x)), 
U*(Z) := Z'(0) + Z(0)(zp/n + /i ), 
:= Z'(tt) + Z(v){ipH x + # ), 

(z, r) = z'(x)y(x) - z(x)Y'(x), 

where Z = [# fc ]*,— is a row vector (t is the sign for the transposition). Then 

(Z, Y) \ x=0 = U*(Z)Y(0) - Z(0)U(Y), (Z, Y) ]x=n = V*(Z)Y(n) - Z{it)V{Y). (15) 
If Y(x, p) and Z(x,p) satisfy the equations £Y(x, p) = and £*Z(x,p) = respectively, then 

^{Z(x,p),Y(x,p)) = 0. (16) 
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Let (p*(x, p), S*(x, p), ip*(x, p) and $*(x, p) be the matrices, satisfying the equation i*Z = and 
the conditions <p*(0,p) = S*'(0,p) = ip*(n,p) = [/*($*) = I, U*(<p*) = S*(0,p) = V*(tp*) = 
= 0. Put M*(p) := $*(0,p). Then 



&(x,p) = S^p) + M*{p)<p*(x,p) = (E/*(V>*))~V*(M- 
According to ([TBI) . ($*(x, p), <3>(x, p)) does not depend on x. Using (fTo^) . we get 

($*(x,p)Mx,P))\x=Q = M (p) ~ M *(P)> ($*(M,$(M)|*=* = 0. 



(17) 



Therefore, M(p) = M*(p). 

Using ffTol) . one can easily show that 



$*'(x,p) p) 
-<p*'(x,p) <p*(x,p) 





"70" 




/ 



Hence, 



(p{x,p) $(x,p) 
<p'(x,p) &(x,p) 



ip{x,p) $(x,p) 
<p'(z,p) &(x,p) 



$*'(x, p) p) 
-<p*'(x,p) <p*(x,p) 



(18) 



We also need the asymptotics for <p*(x,p) and <£>*(x, p). 

Lemma 3*. Suppose x G (0, 7r), z/ = 0, 1, |p| — >• oo; then 
(i) for p G 6^, we have 



<p*Xx,p) 

r iu \x,p) 



-i P y 



2 

{i P y 



exp(-ipx)(I + hi)P*(x) + Oip 1 " 1 exp(\r\x)) 



exp(-tp(TT - x))(I - Hx)P?(x) + 0(p"- 1 exp(|r|(vr - x))) 



&<y\x,p) = -{-ipY~ l exp(ipx){I + h 1 )- 1 {P?{0))- 1 P?{x) + 0{p u - 2 exp{-\r\x)). 
(ii) for p G Qj, we have 

■exp{ipx){I - hJPlix) + 0(p u ~ 1 exp(\r\x)), 

exp(zp(7r - x))(I + ffi)P!*(aO + 0{p v - 1 exp(|r|(7r - x))), 
$* W (z,p) = -(2p)^ 1 exp(-zpx)(/-/n)- 1 (P , *(0))- 1 P , *(x) + 0(p^ 2 exp(-|r|x)). 

Proof of Theorem 1. Consider the problems L and L with the Weyl matrices M(A) = M(A). 
Note that according to asymptotics for M(A) of Lemma [31 we straightway get 



<P* M 


>>p) = 


(v)" 


^*(") 


>>p) = 


Hp) 

2 




>,p) = 


Hp) 



Now we consider the block-matrix V(x,p) = [P jk{ x i p)]j ,k=i,2 defined by 



<p(x,p) $(x,p) 
<p'(x,p) &(x,p) 



P(x,p) 

Taking (TT8"1) into account, we calculate 



<p(x,p) $(x,p) 
<f/(x,p) &(x,p) 



(19) 



(20) 



P n (x, p) = ^ (x, p)$*' (x, p) - (x, p)(p*' (x, p) , 

P j2 (x,p) = ®U- l \x,p)<p*{x,p) -^- l Xx,p)&{x,p). 



(21) 
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Applying Lemmas [3] and 3* and using f lT9|) . for p G 0^~, |p| — >■ oo, x G (0, ir) we calculate 



^P + (a;)( J Pr(0))- 1 Pr(^) + ^ 

1 1^, , " l 



n 2 (x, P ) = — p;(o;)(p;(o))- 1 p;(x) - — P + (x)(P|*(0))- 1 P;*(x) + 0(p~ 2 ). 

Similarly, for p G Qj , \p\ —¥ oo, x G (0, 7r) 

P n (x,p) = Ip.ix^PTm^PTix) + ^'(^(P^O))- 1 ^^) + 0(p" 1 ), 
P 12 (x,p) = -± P -( x )(p-(0))-iP:( x ) + ^-P_(x)(Pr(0))- 1 ^r(a;) + 0(p- 2 ). 

Using Lemma [2], we obtain for p G U 6^, \p\ —> oo, x G (0, 7r), that 

Pn(x,p) = + 0(p- 1 ), P 12 (x,p) = p- x A(x) + 0(p- 2 ), (22) 

where 

fi(z) := ~ (p_(x)P;(x) + P+(x)P!(s 

A(x) : = 1 (p_( x )P;(x) -P + (x)P*(x) 
Substituting (H3J) and (JED into ([21]), we get 

p n = <^*' - S<p*' + <p(M* - M)tp*, 

p 12 = Sfi* - <pS* + <p(M - M*)<p*, 

Note that for each fixed x G (0,7r), the matrix functions <p, <p*, S 1 , 5* and they derivatives 
with respect to x are entire in p of order 1. Since M*(p) = M(p) = M(p), it follows for each 
fixed x G (0, 7r) that the entries of Vn(x, p) and Pi2(x, p) are entire functions in p of order not 
greater than 1. Together with f[2"2"j) this yields Pn(x, p) = f2(x), V\i{x,p) = 0, A(x) = 0. 
Differentiating £7(x) and using (j4j), we get 

Q'(x) = i(A(x)Qi(x) - Q x (x)A(x)) = 

for almost all x in [0, 7r]. Therefore, Q(x) = f2(0) = J, x G [0, 7r]. Thus, Pn(x, p) = J, x G (0, 7r). 
By virtue of (I2"0l) we have <p(x,p) = <p(x,p), $(x,p) = $(x,p) and consequently, L = L. □ 
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